Correlation effects in total energy of transition metals and related properties 
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We present an accurate implementation of total energy calculations into the local density ap- 
proximation plus dynamical mean-field theory (LDA-I-DMFT) method. The electronic structure 
problem is solved through the full potential linear Muffin-Tin Orbital (FP-LMTO) and Korringa- 
Kohn-Rostoker (FP-KKR) methods with a perturbative solver for the effective impurity suitable for 
moderately correlated systems. We have tested the method in detail for the case of Ni and investi- 
gated the sensitivity of the results to the computational scheme and to the complete self-consistency. 
It is demonstrated that the LDA-I-DMFT method can resolve a long-standing controversy between 
the LDA/GGA density functional approach and experiment for equilibrium lattice constant and 
bulk modulus of Mn. 

PACS numbers: 71.15.Nc, 71.20.Be, 71.27.+a 



I. INTRODUCTION 

The state-of-the-art technique for calculating the 
electronic structure of materials is density functional 
theorji^ in its local density approximation (LDA). 
However, despite numerous impressive successes, it 
faces serious difficulties for strongly correlated systems 
such as Mott insulators, heavy fermion systems, high- 
temperature superconductors, itinerant electron mag- 
nets, and many others. Some of these difficulties were 
overcame by merging LDA-based first-principles elec- 
tronic structure calculations with the dynamical mean- 
field theory (the LDA-I-DMFT approach'^-'*; for review 
see Refs. USQIl). Most of the works done by this 
method deal with spectral properties of strongly corre- 
lated systems. At the same time, correlation effects are 
sometimes of crucial importance to describe also cohe- 
sive energy, equilibrium lattice constant and bulk modu- 
lus, as demonstrated for the cases of plutonium^iiiS and 
ceriumiiii^. After these first attempts, the need of more 
systematic implementations and investigation of the nu- 
merical problems related to total energy evaluation in the 
LDA-I-DMFT scheme arose. Recently Pourovskii et ali^ 
have presented an interesting comparison between the 
correlation effects in the basic DMFT cycle (convergence 
in the local self-energy) and in the fully self-consistent 
DMFT cycle (convergence in the local self-energy and in 
the electron density). Two test-cases have been studied 
with this new implementation: the 7 phase of metallic 
cerium and the Mott insulator Ce203. Both of them arc 
close-packed /-electron systems and they can be studied 
through the atomic sphere approximation within the lin- 
ear muffin-tin orbital method (ASA-LMTO) and through 
the Hubbard-I solver—. 

Up to now all the LDA-I-DMFT studies of the ground 
state properties of strongly correlated systems concerned 



materials with rather localized /-electrons. Here, we 
are interested in materials where the correlation effects 
are less dramatic and where the failures of the density- 
functional theory pertain only some specific proper- 
ties. The late transition metals Mn, Fe, Co, and 
Ni are correlated systems, and the LDA-I-DMFT ap- 
proach was successfully applied to describe their spec- 
tral propertie o-^^i^^'^^i-^'^'^^i^^'^°i^-^i^^ as well as their mag- 
netic propertiesi^S. In particular, the DMFT was 
implemented into full-potential Korringa-Kohn-Rostoker 
method (FP-KKR)i^ and full-potential linear muffin-tin 
orbital method (FP-LMTO )^^ to allow corresponding 
studies. 

In the present paper we extend the previous imple- 
mentations to calculate the total energy of the electronic 
system whitin the LDA-I-DMFT scheme. While there 
exists already another LDA-I-DMFT code based on FP- 
LMTO and able to calculate ground state properties of 
strongly correlated materialsi^, we must emphasize that 
it is the first time that this happens for an LDA-I-DMFT 
code based on FP-KKR. This is particularly appealing 
to analyze disordered alloys systems, for which FP-KKR 
in combination with coherent potential approximation 
(CPA) alloy-theory is well-known to be particularly ef- 
ficient and reliable. 

Here we use these implementations to study the to- 
tal energy and related properties of 3c? transition met- 
als. First, we present computational results for Ni which 
plays the role of "drosophila fiy" for the LDA-I-DMFT 
method and where the most detailed comparison of the 
theory with experiment was done^5,iii. After calculations 
of photoemission, optical and magnetooptical spectra, 
magnetization, magnetic susceptibility, orbital magnetic 
moments, bulk and surface spectral densities (see previ- 
ous works cited above), the present calculation of cohe- 
sive energy, equilibrium lattice constant and bulk mod- 
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ulus completes its basic physical description within the 
LDA+DMFT approach. 

Comparing the results of the full-potential KKR and 
LMTO calculations we address the question about sensi- 
tivity of the LDA+DMFT description to the band struc- 
ture method used. This is nontrivial since different meth- 
ods use different basis sets which are truncated in any real 
calculations. We have found that actually the results are 
very close which support its reliability. While correla- 
tion effects in ground state properties of Ni are quite 
small, they are accurately described within our scheme 
which confirms the usefulness of the LDA-I-DMFT for 
not only strongly correlated but also for moderately cor- 
related systems. We have checked also the importance 
of the full charge self-consistency and found that in the 
case if Ni these effects are not very essential. 

Then, we consider the case of Mn where, among 
all transition metals, the largest discrepancy between 
the LDA or GGA predictions for the lattice constant 
and bulk modulus and the experimental data takes 
place^ii^i^ which is considered to be an indication of 
the strongest correlation effects among 3d metals^i^l. 
We show that the LDA+DMFT method does allow us to 
solve this problem and to describe in a very satisfactory 
way the energetics of Mn. 

II. FORMULATION OF THE PROBLEM 

All the standard approaches for calculating the elec- 
tronic structure of strongly correlated materials are based 
on the choice of a set of orbitals that are described not 
accurately enough in the standard DFT-LDA method 
which is supposed to be improved. We call them "corre- 
lated orbitals" and indicate with |R, <^), where R is the 
vector specifying the Bravais lattice site and the ^ is an 
index that enumerates the orbitals within the unit cell of 
the crystal. The choice of {|R, £,)} is dictated by physical 
motivations for the problem under consideration and al- 
ways implies some degree of arbitrariness (see the discus- 
sion below). Usually the correlated orbitals are derived 
from d or / atomic states and the index ^ stands for the 
atomic quantum numbers /, m, a. Natural choices can be 
Linear Muffin-Tin OrbitalsiS, or Wannier functions^i^. 
Apart from the atomic states, hybridized orbitals can 
also be chosen depending on the problem. For example 
in the transition metal oxides the crystal field splits the 
LDA bands in two distinct groups, well separated in en- 
ergy and suitable to be determined through downfolding 
of the original problem via the NMTO approach'^". 

After having decided the set {|R, ^)}, we correct 
the standard DFT-LDA Hamiltonian with an additional 
Hubbard interaction term^ that explicitly describes the 
local Coulomb repulsion U for the orbitals in the set: 

(1) 



This is the so-called LDA+U Hamiltonian and an im- 
portant remark has to be made concerning the mean- 
ing of the matrix elements C^^i,^2,C3,C4- We should not 
think of them as generic matrix elements of the bare 
Coulomb repulsion, but more as the matrix elements 
of an effective interaction introduced to give the cor- 
rect description of the low energy excitations (to de- 
scribe broader energy scales the U term should be, in 
general, energy dependent"^^). In these terms we have 
to consider the LDA+U Hamiltonian as being derived 
from completely ab initio density functional approach. 
While in principle this is possible through many methods, 
e.g. constrained density functional theor y'^^i'^^ , extrac- 
tion from GW results^l^ it is a common practice to eval- 
uate the matrix elements [7^^^2.43.44 using semi-empirical 
procedures''*''^'*. This may seem to be inadequate, since 
the strength of the effective Coulomb interaction should 
depend on the set of correlated orbitals, being strictly 
connected to a mapping of the original electronic Hamil- 
tonian into the Eq. ([T]); however if the orbitals {|R, 1^)} 
are chosen appropriately, the results are quite stable with 
respect to this ambiguity, as it was first noticed for the 
LDA+U method^^. 

If the correlated orbitals are atomic-like ones (with the 
quantum numbers Z, m, cr), we can express^ the Coulomb 
parameters in terms of Slater integrals i^" like 
21 

^Cl,?2,43,«4 =^^an{£,i,£,3,^2,£,4)F"', (2) 
with the coefficients a„ defined as 

an(Ci,6,6,C4) = ^^ E (aii^«,i6)(6ii;:jC4), 

q— — n 

(3) 

where the terms (Cil^gl^a) and (C2|^g|^4) are integrals 
over products of three spherical harmonics. In the rest of 
the paper we will be interested in 3d electrons, therefore 
we can limit our discussion to these. For d electrons 
there are only three Slater parameters (F°, and F*) 
and they can be easily linked to the Coulomb parameter 
U and the Stoner parameter J as2^ 



The ratio F^/F"^ is assumed to correspond to the atomic 
value and for 3d electrons it is approximately equal to 
0.625. In the rest of the paper we will use the two real 
values U and J to specify the Coulomb matrix elements. 

III. DYNAMICAL MEAN-FIELD THEORY 

The LDA+U Hamiltonian defines an "effective" Hub- 
bard model, and its solution represents a complicated 
many-body problem. The strategy of the spectral den- 
sity functional theory^ is the same of DFT or Baym- 
Kadanoff theory (or more generally of every Weiss mean- 
field theory): we specify a main observable quantity and 
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we map the original system into a system with less de- 
grees of freedom under the condition of conserving the 
expectation value of the main observable. In DFT and 
Baym-Kadanoff theory the main observables are respec- 
tively the total electron density p(r) and the one electron 
Green's function G(z), namely 



LOCAL PROBLEM 



G(z) 



t{z) 



(5) 



where z is the energy in the complex plane, jj, is the chem- 
ical potential, Hlda plays the role of the unperturbed 
Hamiltonian ("hopping"), and '^{z) is the self-energy op- 
erator, which in many-body theory reproduces the effects 
of the interactions. In spectral density functional theory 
the main observable is the local Green's function at the 
site R: 



Gr(z) = PrG(z)Fr, 



where 



Pi 



R 



i 



(6) 



(7) 



is the projection operator to the correlated subspace be- 
longing to site R. 

As in density functional theory, where we make ap- 
proximations as LDA or GGA, in the framework of the 
spectral density functional theory the corresponding ap- 
proximation is the dynamical mean-field theory. In the 
DMFT the self-energy is considered to be purely local. 
In terms of matrix elements on the correlated orbitals at 
the two sites Ri and R2, this means that 

(Ri,a|sW|R2,C2) = -Jr^r, (6|Sr, • (8) 

The assumption of a purely local self-energy Sr(2:) al- 
lows us to focus only on the single lattice site R. Con- 
sequently we can replace all the other sites of the lat- 
tice with a self-consistent electronic bath (or "dynamical 
mean- field") C^q"^(R, z), whose role is analogous to the 
Weiss mean-field used in statistical mechanics. What 
we have obtained is a problem of an atomic site em- 
bedded into the fermionic bath: in many-body physics 
this system is known as multi-band Anderson Impurity 
Model. While we do not have obtained an Hamiltonian 
that describes the mapping system, we can easily write 
down the effective action 5, so that the problem is fully 
determined. All the mathematical details, the explicit 
formulas and a more detailed description of the DMFT 
equations can be found in Refs. @ or 0. 

The Anderson Impurity Model has been widely studied 
in the many-body literature and its solution can be ob- 
tained through many different techniques, usually named 
"solvers" in the DMFT community. At the present time 
no solver has succeeded to become the "standard ap- 
proach" of the LDA-I-DMFT scheme, but the technique 
to be used is every time decided with respect to the 
strength of the correlations and the degree of accuracy 
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FIG. 1: (Color on-line) Schematic representation of the iter- 
ative procedure to follow in the LDA-I-DMFT scheme. As a 
first step the DFT-LDA problem is solved and a ground state 
electron density p(r) is obtained. From p(r) we can extract 
the matrix elements of the single-particle LDA Hamiltonian, 
and then build the one-electron Green's function G(k, ici;„) 
at the Matsubara frequencies iujn- Now the basic DMFT cy- 
cle starts: the Green's function G(k, itj„) is projected onto 
the correlated orbitals, defining the bath Green's function 
C/(7^(R, iijJn) of the Anderson Impurity Model by means of the 
Eq. ([SJ . The solution of the local problem through one of the 
available "solvers" leads to a self-energy function ER(iaj„). 
After a back projection to the LDA basis set, a new one elec- 
tron Green's function G(k, ia;„) and a new chemical poten- 
tial \x are calculated. The procedure is repeated iteratively 
until convergence in the self-energy and the chemical poten- 
tial. Once the convergence of the basic DMFT cycle has been 
reached, a new electron density p(r) can be calculated from 
G(k, icij„). This is the fully self-consistent cycle and should 
be continued until convergence in p(r). 



desired. In the case of systematic simulations, as for 
example the total energy calculations reported in this 
paper, another important factor to consider is the nu- 
merical efficiency. In both the implementations of the 
LDA+DMFT scheme discussed here the spin-orbit spin- 
polarized T-matrix ffuctuation-exchange (so SPTF— i^) 
solver has been used, being reliable and efficient for mod- 
erate strength of the correlations (J7 < IF/2 where W is 
the bandwidth of the localized orbitals). 

Once the effective impurity problem has been solved 
and a self-energy Sr(z) has been obtained, there is an 
apparent change of the number of particles, so that the 
chemical potential /i has to be updated. Furthermore 
a new electronic bath ^(^^(R, z) is defined through the 
inverse Dyson equation: 



^o-1(R,z) = Gr1(z) + Er(z). 



(9) 



Now we can iterate the procedure described above until 
convergence of the self-energy and the number of par- 
ticles (or chemical potential). This is the basic DMFT 
cycle and is schematically reproduced in Figure 1. From 
the same Figure, we can also notice that, if the correla- 
tions are strong, the differences in the population of the 
Kohn-Sham orbitals lead to a new electron-density /o(r). 
In this "full self-consistent cycle" also the convergence 
of p(r) has to be reached. In the present paper two im- 
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plementations of the LDA+DMFT scheme are used: the 
first one works only within the basic DMFT cycle, while 
the other one uses the fully self-consistent cycle. There 
are mainly two reasons why we have compared these two 
different implementations. First of all we want to study 
3d states: they are not extremely localized and conse- 
quently the effective Coulomb repulsions between them 
are not very strong. Then we can reasonably suppose 
that the changes in the electron-density are small, and 
a measure of this is given by the comparison of the two 
different codes. In second place one of the aims of our 
study is to investigate which numerical precision can be 
obtained for the LDA-I-DMFT scheme, also in compar- 
ison to the modern DFT packages, so to allow reliable 
calculation of sensitive quantities as equilibrium atomic 
volume and bulk modulus. With respect to this issue, we 
should stress that implementing full self-consistency over 
the charge density is a very delicate task that can bring 
additional numerical errors. Moreover the computational 
effort can rise considerably, limiting the applications of 
the LDA-I-DMFT scheme to systems with a few atoms 
per unit cell. 

Before presenting the total-energy functional adopted 
in the LDA-I-DMFT scheme, a final remark has to be 
made. Since the LDA-I-U Hamiltonian is constructed 
with an additional term that is already contained in 
the original electronic Hamiltonian, we have to remove 
from the self-energy those contributions already calcu- 
lated in the LDA. Unfortunately, there is no way to 
establish exactly a correspondence between approxima- 
tions within the density functional and the Green's func- 
tion (Baym-Kadanoff ) functional, then we have simply to 
"guess" which diagrammatic contributions are included 
and which ones are not. For treating metals the most 
common choice of the "double counting" correction is the 
static part of the self-energ y "'^^i^^ . In the present paper 
we adopt the double counting of Ref. i.e. we treat the 
static contribution to the self-energy as in the LDA-I-U 
method with around mean-field (AMF) double counting, 
while the other contributions to the self energies become 

S^,,5,(z) = %,«,(z)-%,e,(S(0)), (10) 

where the average has to be determined over the orbital 
indices separately per spin channel. This choice is due to 
the fact that the LDA exchange-correlation potential is 
an orbitally averaged quantity and has proven to be very 
succesfuU in describing the transition metals. 

IV. TOTAL ENERGY FUNCTIONAL 

In the previous section we have presented the equa- 
tions that define the LDA-I-DMFT scheme in terms of 
local problem and self-consistent bath. These equations 
can be obtained with many different techniques^, but 
in perspective of total energy calculations we have al- 
ready adopted the point of view of the spectral density- 
functional theory of Savrasov and Kotliar. In a series of 



paper^»i2. they have introduced a functional of both the 
total electron density p(r) and the local Green's func- 
tion Gr(z) for the correlated orbitals. It is important to 
emphasize that these quantities are independent, in the 
sense that they cannot be reconstructed from each other. 
Furthermore notice that, in this framework, the arbitrari- 
ness of the basis set of the correlated orbitals is contained 
in G'r(z). Following standard methods of quantum-field 
theory the functional is constructed introducing source 
terms for p(r) and Gr(2); then the variational procedure 
is applied to the functional with respect to these sources. 
Without presenting the mathematical details (see refer- 
ences above), we obtain the following expression for the 
zero-temperature limit of the total energy: 

E = Elba [p(r)] - £k. + TtIHldaG] + {Hu) (H) 

where Hjj indicates the two-particle term in the LDA-fU 
Hamiltonian ([T|), and the primed sum is over the occu- 
pied states. Here and in the following the symbol Tr 
indicates the one-electron trace for a generic represen- 
tation and the sum over the Matsubara frequencies iuj 
for finite temperature many-body formalism. We assume 
that the temperature effects can be taken into account 
only via summation over the Matsubara frequencies and 
in the DFT part only weak temperature dependence via 
the Fermi distribution function is taken into account^. 
This corresponds to neglect the temperature dependence 
of the exchange-correlation potential and it is a standard 
procedure in electronic structure calculations of real ma- 
terials. These effects are irrelevant for the cases under 
consideration where the main temperature dependence 
is due to spin fluctuation&ii. 

We notice that the total energy within the 
LDA-I-DMFT scheme is not simply the expectation value 
of this Hamiltonian, but it consists of several terms, 
in analogy to the expressions of the usual DFT. The 
first term E^da [p(r)] contains four different contribu- 
tions, namely the ones due to the external potential, 
the Hartree potential, the exchange-correlation potential 
and the sum of the Kohn-Sham eigenvalues. However 
in the spectral density functional theory the Kohn-Sham 
eigenvalues should be re-occupied with respect to the de- 
scription given by the total Green's function. Then we 
should remove the bare Kohn-Sham eigenvalues sum (sec- 
ond term of Eq. pTjl ) and substitute it with Tr [H^daG] 
(third term). Moreover notice that Eld A [p(r)] depends 
only on the total electron density, so it does not need to 
be recalculated if the basic DMFT cycle is applied. In 
the case of the fully self-consistent cycle, the calculation 
is straightforward, since it uses the same LDA-DFT ma- 
chinery. This point will be analyzed in more details in 
the section concerning the FP-KKR implementation. 

Finally we can evaluate {Hjj) through the so-called 
Galitskii-Migdal formul a^^i^° , an elegant way to rewrite 
the expectation value of a two-particle operator in terms 
of a one-particle operator as the Green's function. This 
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formula is based on the fact that for an Hamiltonian 
H — Hq + Hif, i.e. the same form of the Hamiltonian 
^ , the equation of motion of the Green's function states 
that 



(12) 



where r is the imaginary time for the finite temperature 
formalism (the formulation for real times and T = is 
completely equivalent). Using the Fourier transform with 
respect to r, we can move to the energy domain 



|:G(r)) = Tr[c^G(^)]. 



(13) 



Furthermore from the definition of the Green's function 
[wl — Hq — E(aj)]G'(a;) = 1, we can rewrite the expres- 
sion above in terms of more manageable operators 

Tr [ujG{uj)] ^ Tr [t{uj)G{uj)] + Tr [HaG{uj)]. (14) 

Then the Galitskii-Migdal formula becomes 



{Hu) = lTT[tG] 



IMPLEMENTATION IN FP-LMTO 



(15) 



We have implemented the total-energy algorithm of the 
previous section in the recently developed LDA+DMFT 
code^, based on the full-potential linear mufhn-tin or- 
bital (FP-LMTO) method code presented in Ref 41 and 
well-known to give accurate description of many solids 
within LDA. The full-potential character of the program 
makes it very attractive for open structures and surfaces, 
and in fact the first applications of our code were focused 
on these systems^^ . Furthermore the use of a small num- 
ber of basis functions as used within the LMTO method 
is particularly efhcient for calculating the Green's func- 
tions, since they require inversions of a matrix in the 
LDA basis set for each Matsubara frequency and k point. 
While we do not want to give a complete survey of the 
equations involved in the FP-LMTO code, the study of 
the total energy problem forced us to develop a more so- 
phisticated method to calculate the number of electrons 
for the given chemical potential. In this implementation 
two basis sets are used: the already mentioned set of the 
correlated orbitals {|R, ^)} and the set of the LDA basis 
functions {|k, x)}- The steps of the LDA-I-DMFT scheme 
that imply moving from the local problem to the lattice 
problem require transformations back and forth between 
these two basis sets. Furthermore we should mention 
that the set {|k, x)} is neither normahzed nor orthogo- 
nal and then the basic algebraic operations involve an 
overlap matrix 



and its inverse 5'""'^, since the dual basis set of {|k, x)} 
does not coincide with the set itself. The number of elec- 
trons in the lattice problem is calculated with the LDA 
basis set and becomes 

= T^E E E ^(k)x.,xiG(k, zo;„)x„^, , (17) 

iuj„ k XI, X2 



where 



G(k, icuri 



{k,Xl\G{iL0n)\k,X2). (18) 



The sum over the Matsubara poles should include infinite 
negative and positive frequencies, but obviously in a com- 
putational scheme the number of frequencies can only be 
finite and then a cut-off value uJmax needs to be chosen. 
Unfortunately, as it is clear from the definition ([5]), the 
Green's functions decay slowly with the energy and then 
a reliable determination of the number of particles would 
require a huge cut-off. There are two ways to reach this 
cut-off: increasing the number of Matsubara frequencies 
or increasing the spacing between them, proportional to 
the temperature T. None of them is a good solution. The 
former would imply too big numerical effort (there is an 
inversion of a matrix with the size of the LDA basis set 
for every Matsubara frequency and every k point), while 
the latter would lead us too far from the T = ground- 
state. In Ref. [13 the problem of the long-decaying tails 
of the Green's function was solved in a rather rudimen- 
tal way, given that the paper was focused on the spec- 
tral properties, which are not as sensitive as the ground 
state properties to the numerical details. In the present 
paper, conversely, we follow the elegant procedure used 
in the LDA+DMFT calculations'^' and adapted to our 
non-orthonormal basis set. The idea is to decompose the 
calculated Green's function as 

G(k, zc.„)^, = G(k, zc^„);jr."x. + G(k, ^cOnrxlx2 > (19) 

where G(k, j(x>„)J" is a analytical function that we 
chose to fit the high-frequency behavior of G(k, icj„)^j ■ 



E 



[G(k, 



"'"■JXl,X2 



G(k,iLL!n] 



an 

Xl,X2j 



= 0. (20) 



On the other hand the numerical part is defined as the 
difference between the calculated function and the ana- 
lytical function 



G(k,z^„)— 



G(k,zc.„)x,,^, -G(k,zc.„)«;',^^, (21) 
has been chosen wisely, is negligible 



and, if G(k, i'^n)xi,x2 

for UJn > iOmax- 

The new problem is to determine G"". Starting from 
the definition (ISl) , we can rewrite the matrix element (1181) 



G(k,iw„)xi.X2 = ( k,x 



5(k), 



(k,Xi|k,X2) 



(16) 



k, X2 



(22) 
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where we have defined the new operator 

Ik(it^n) = /-ii - fiLDA - S(iu;„) . (23) 

Let's consider S(ia;„) — corresponding to the first it- 
eration of the LDA+DMFT cycle. In this case the oper- 
ator pS)) does not depend on the Matsubara frequencies 
and is Hermitian; consequently it has real eigenvalues 
A}^ and the eigenvectors \X^) can be chosen to form an 
orthonormal set. By expanding A]s^(iLUn) in its spectral 
representation, the Eq. (|19p becomes 



^-^ iujn — All 



We have finally reduced the original sum to two terms 
that we can calculate with high precision. The numerical 
term is simply calculated as a sum for positive frequencies 
up to Umax ■ The sum for negative frequencies is obtained 
using the symmetry of the Green's function 



G(k,-iujn)xux2 = [G'(-k,it^„)x2.xi 



(25) 



The analytical term can be summed through standard 
many-body techniques: 



EE 



(k,xi|Xk)(Xk|k,X2) 



Ajl 



(k,Xl|X^J(Xk|k,X2) 



1 



(26) 



In comparison with Ref . |l3| finding eigenvalues and eigen- 
vectors of Ak is slightly more cumbersome here: due to 
the non-orthonormality of the basis set we have to solve 
a generalized eigenvalue problem. However using the 
fact that the overlap matrix is positive definite, through 
Cholesky decomposition^^ of S the problem can be re- 
duced to a usual eigenvalue problem through a few alge- 
braic operations. 

When the DMFT self-energy assumes finite values, the 
operator Ak(i(^n) is different at every Matsubara fre- 
quency, and then we need to use some approximation. 
Luckily in many-body theory the analytical properties of 
the self-energy operator are the same as for the Green's 
function. Therefore we can assume the following asymp- 
totic behavior for high frequencies-'^'^: 



Y^asyn 



lUJ 



(27) 



where E^*''* and S^^J"" are obtained from the real and 
imaginary part of S at the cut-off value oJmax- While a 
higher oJmax will always give a better fit, the real part 
of the self-energy converges to E'**"* as and then 

we do not need a very high cut-off. Furthermore for our 
purposes of evaluating the frequency sum in Eq. (|17p , we 



can keep only the dominant term E'**"* , and E"''!'™ turns 
to be unimportant. Again the operator (|23p becomes 
Hermitian and independent on the Matsubara frequen- 
cies, so that the same procedure described above can be 
applied. 

The implementation of this algorithm in the FP- 
LMTO code improved the precision in the determination 
of the number of particles by about two orders of magni- 
tude in the worst cases (corresponding to a large number 
of LDA basis functions that increases the numerical er- 
ror on the eigenvectors). The method used in Ref. [13 
was rather similar to the one presented above, but had 
a much simpler implementation. Instead of considering 
the asymptotic behavior of every Green's function in Eq. 
()17|1 . the sum over the intermediate indices XI1X2 a-nd k 
was done, and then the asymptotic behavior of the result- 
ing function was considered. While this approximation 
can appear too crude, the precision on the number of par- 
ticles is about 10^"^ particles for every electron involved in 
the problem. On the other hand it was computationally 
very efficient, since the generalized eigenvalue problem 
was reduced to the determination of a pure number. 

After having improved the precision in the determina- 
tion of the number of particles, we can pass to the imple- 
mentation of the total energy formula ([TT|) . As we have 
already seen the first two terms can be obtained from 
the density-functional part of the LDA-f DMFT scheme. 
The third term, corresponding to the reoccupation of the 
Kohn-Sham orbitals, requires again the evaluation of a 
sum over all the Matsubara frequencies 

Tt[HlDaG]=T^'Y^ E HLDA{^)x2,XiG{^ii^n)xuX2- 
iuJn k Xl>X2 

(28) 

Besides the presence of different matrix elements, Eq. 
(p8)) is completely analogous to Eq. (fT7|) . therefore the 
sum can be done by applying the same procedure used 
above. Finally we have to evaluate the Galitskii-Migdal 
contribution {Hu}- Given that in the LDA-hDMFT 
scheme the self-energy is local, the trace in Eq. (flS)) 
can be restricted to the correlated orbitals. Further- 
more, using the fact that in the SPTF solver we work 
with quantities in both the frequency and (imaginary) 
time domains, we can express the trace in terms of the 
complex Fourier transforms. For this purpose, it is most 
convenient to separate the static and the dynamic parts 
of the self-energy. Analogously to Eq. ((27|) . we have 

(29) 



However now no fitting is necessary: once S''*'^' is de- 
termined, Ti{iuinY^^ contains all the differences with the 
calculated function S(z<jj„). We can then write 

= ^^E E +E(zc.„)^fl]G(zc.„)6,6. 

(30) 

The first term at right hand side can be easily Fourier 
transformed and reduced in terms of occupations of the 
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local orbitals 



(31) 



the second term requires the evaluation of the Fourier 
transform of a product, leading to a convolution. In sum- 
mary we can express Eq. (j30p as 



- f dTnT)f-^G{~T\,,^,]. (32) 
Jo 



DMFT 
(FLEX) 




VI. IMPLEMENTATION IN FP-KKR 

The same total energy algorithm of the previous sec- 
tions was implemented in the FP-KKR code described 
in Ref. [H, being an extension to the full-potential case 
of the earlier ASA implementationi^. Besides the ad- 
vantage of being one of the very few fully-self consistent 
implementations of the LDA-I-DMFT scheme, the for- 
malism on which the FP-KKR code relies makes it par- 
ticularly attractive to study complex problems as orbital 
polarizations^, photoemission spectroscopy through the 
one-step model^i, or disordered alloys systems through 
CPA43. drawback to the flexibility of LDA+DMFT, 
in FP-KKR we have a high computational cost that can 
make it inconvenient to perform extensive simulations, 
e.g. determination of total energy curves as functions of 
the crystal parameters, compared to other simpler meth- 
ods. 

Without presenting a complete survey of the equations 
involved in the FP-KKR method, we should mention that 
it is based on the multiple scattering theory which al- 
lows to decompose the total one-particle Green's func- 
tion into the single scattering matrix t^{e) which con- 
tains the information about each single scatterer R (i.e. 
atomic site) and the structure constants matrix (e) 
which contains the information about the geometrical ar- 
rangement of the scatterers in a solid. All the ingredients 
are calculated in the basis of the four-component energy- 
dependent regular and irregular solutions of the rela- 
tivistic Kohn-Sham-Dirac equation s — '™ . Corresponding 
DMFT self-energy is included into corresponding quasi- 
particle Dirac equation, e.g. in contrast to the other 
LMTO based LDA-hDMFT implementations the corre- 
lation effects are included directly into the single site t 
matrix as well as into the wave functions at the same 
time. 

In order to construct the bath Green's function needed 
as an input for the DMFT solver, the localized Green's 
function is calculated by projecting the total Green's 
function onto the correlated atomic site. The multiple 
scattering formalism provides the natural choice of the 
projectors which are nothing else as the regular single- 
site solutions of the Kohn-Sham-Dirac equations. The 
projection functions are taken at the fixed energy, which 



FIG. 2: (Color on-line) Illustration of the KKR-(-DMFT 
scheme: blue semicircle is the complex energy path used 
by KKR to calculate the Green's function. After the bath 
Green's function G is obtained, it is analytically continued 
onto the imaginary axis (red) to calculate the self-energy via 
the SPTF impurity solver. The latter is analytically extrap- 
olated back to the semicircle. 

corresponds to the center of mass of the band and is re- 
calculated at each iteration. 

Different from FP-LMTO, the FP-KKR code works 
with the Green's functions from the beginning, so that 
the merging between the LDA part and the DMFT part 
in the LDA-I-DMFT scheme does not require a change 
of representation of the electronic states. In practice, 
as shown in Figure [2] the LDA Green's function is eval- 
uated on a semicircular contour in the complex plane, 
but the SPTF solver works with the Green's functions 
on the Matsubara frequencies. The analytical continu- 
ation of the self-energy from the Matsubara axis to the 
semicircular contour is done through the Fade approxi- 
mants method, and this could introduce small numerical 
errors. While no problem was observed in all the previ- 
ous studies for which this code has been used, we have 
considered that the determination of the energetic land- 
scape requires more precision than spectral properties. 
For that reason we have checked this point carefully. As 
expected we have found a small numerical noise, but in 
practice its effects on the ground state properties of the 
transion metals studied here were negligible. 

In FP-KKR the total energy functional (fTT|) can be 
rewritten in a slightly different form. In the LDA contri- 
bution to the LDA-I-DMFT total energy we can explici- 
tate the standard terms that are going to be summed in 
density functional theory. We have 

Elba W)] ^ J ^ de eNR{e) ~ J dr p{r)VR{r)+ 

(33) 

where i? is a given region of the space defined in the 
unit cell and Nr is the number of electrons in the space 
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R. This occupation number is obtained from the un- 
perturbed LDA Green's functions, that is from the in- 
formation specified by p{r). This is the reason why the 
functional dependence of the energy above is restricted 
to the only electron density. If we calculate the num- 
ber of electrons with the full DMFT Green's func- 
tions, we obtain that the first term of the functional 
the so-called "band energy", becomes exactly the 
term at left hand side of Eq. (fT4)) . Then, renaming 
the quantity f33|) calculated with this new occupation 
as Eld A [p(r), it is straightforward to rewrite the 

functional pT|) as 

Elda+dmft = Elba [p{r),G{Lo)] - {Hu) (34) 

where the Galitskii-Migdal term has to be subtracted 
since it is already accounted for twice within the band 
energy. 

The evaluation of the formula requires only the 
calculation of the Galitskii-Migdal energy ([15]), since the 
band energy results from the DFT part of the FP-KKR 
code. While it could be simpler to evaluate the Galitskii- 
Migdal correction directly on the local problem through 
the formula (j32p . we prefer to work again on the semicir- 
cular complex contour, retaining to the same formalism 
for both the contributions to the total energy. Then we 
calculate 

^^^^ = "^^"'^ / (35) 

The integration is performed over the contour starting 
close to the real energy axis at the bottom of a valence 
band and ending at the Fermi energy. It turned out that 
this procedure is numerically more stable than evaluation 
of Galitskii-Migdal correction using integration over the 
Matsubara frequencies. 

VII. FCC NI 

Bulk fee Ni is a sort of standard test-case for every 
approach to strongly correlated materials. For this rea- 
son it has been chosen as first application for the imple- 
mentations presented above. The interest of the DMFT 
community in Ni startedi^ with the explanation of the 
famous 6 eV satellite observed in photoemission experi- 
ments, but missing in all DFT calculations. Afterwards 
spectral properties of bulk Ni were studied through differ- 
ent LDA+DMFT implementation o-^^'^^ and also through 
the GW-I-DMFT calculations^. All these studies con- 
firmed the correlated nature of the Ni satellite. Fur- 
thermore recent LDA-I-DMFT based calculation of the 
one-step model photoemission spectrum showed a very 
good quantitative agreement with experimental photoe- 
mission data^^. Along with these spectral features, the 
LDA-I-DMFT method has been applied to the finite- 
temperature magnetism'^^ of Ni, showing the existence of 
local moments (unordered above the Curie temperature). 



i.e. another clear sign of strong correlation. Nevertheless 
we should consider that the DFT scheme is not focused 
on the excitation spectrum, but mainly on the electron 
density. Given the the LDA-I-DMFT scheme and the 
Hamiltonian ([T]) are explicitly build for the correct de- 
scription of the low-energy excitations, it appears natu- 
ral that this scheme performs convincingly better than 
simple density functional theory. Conversely DFT gives 
a reasonable description of all ground state properties 
of Ni and the agreement with the experimental data be- 
comes almost perfect if GGA is used^i^i^. Moreover, 
in contrast with the other late transition metals, the in- 
clusion of the spin polarization in the calculations for fee 
Ni is not strictly necessary, surely due to the small mag- 
netic moment (/i ~ 0.6) acquired^ at the equilibrium 
structure. Finally, a recent accurate study of the orbital 
and spin polarization of the late transition metals^'^ em- 
phasized that the DMFT corrections to the DFT-LDA 
values for Ni are really minor, while still improving the 
description of the material. 

With reference to the previous arguments, it appears 
necessary to clarify the reasons behind our interest in the 
energetics of fee Ni, where the correlation effects are ex- 
pected to have a moderate role. First of all it is important 
to complete the picture outlined above: excitation spec- 
trum, magnetism, photoemission spectrum, surfaces, or- 
bital polarization and now ground state properties. This 
study can help in understanding how correlated fee Ni 
is^ and which deficiencies of the DFT-LDA technique are 
due to a single-particle approximation of the exchange- 
correlation potential and which ones are due to the intrin- 
sic meaning of the Kohn-Sham quasiparticles as fictitious 
excitations. In second place Ni represents a good test- 
case to prove the ability of the LDA-f DMFT scheme to 
catch moderate correlation effects in a real material. In 
fact we know that the LDA-I-DMFT scheme relies mainly 
on two different approximations: finite number of near- 
est neighbors (due to the locality of the self-energy) and 
non-exact solver. Therefore it is interesting to check how 
dominant are the errors connected to these approxima- 
tions for effects that are expected to be rather small. Fur- 
thermore a third important question concerns the role of 
the full self-consistency in the DMFT cycle. Previous 
studiesiS for Ce2 03 and 7-Ce have shown, quite surpris- 
ingly, small differences between the ground state proper- 
ties for the basic and fully self-consistent DMFT cycles. 
Given that these systems involve valence electrons much 
more localized than the ones of Ni, in our case we expect 
negligible differences, at least in the range of "accept- 
able" Hubbard U. This would represent a further vali- 
dation of our previous studies^^ of bulk and surface Ni, 
founded on the basic LDA-I-DMFT cycle, limiting the 
necessity of the full cycle to the most sensitive quanti- 
ties like photoemission spectrum^! and magneto-optical 
properties^. Finally, a last question investigated for fee 
Ni concerns the compatibility between different imple- 
mentations: can different codes with different choices of 
the correlated orbitals give comparable results? 
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FIG. 3: (Color on-line) Energy vs lattice constant curves 
for fee Ni in the DFT-LDA scheme and in the LDA+DMFT 
scheme based on the FP-LMTO (top) and FP-KKR method 
(bottom). The zero of the energy of each curve is set to its 
own minimum value Eo and three chosen values of U are pre- 
sented (T — 400 K). The experimental value of the lattice 
constant is indicated by the arrow. 



To investigate all the various points outlined in the 
previous paragraph, we performed LDA-f DMFT simula- 
tions of fee Ni for various lattice constants starting from 
a = 6.2 a.u. and up to a = 7.4 a.u.. We treated 3d, 4s 
and 4p electrons as valence electrons. For the FP-LMTO 
simulations, the description of the valence electrons in the 
interstitial space between the muffin tin spheres requires 
LMT-Orbitals with different tail energies, whose number 
depends on the degree of localization-delocalization of the 
electrons: three tails were used for 4s and Ap electrons, 
only two tails for 3d electrons. The set of the corre- 
lated orbitals was build from the LMT-Orbitals, consid- 
ering only the part contained into the muffin-tin sphere 
at a given linearization energy, the so-called "head of 
the LMTO"—. Convergence on the total energy with re- 
spect to the k-mesh lead to a minimum number of 4913 
k-points used in the three dimensional Brillouin zone. A 
simulation has been considered converged if the energy 



difference for two consecutive iterations has been at least 
smaller than 0.1 meV. As far as possible same settings 
were used for the FP-KKR simulations with the excep- 
tion of the set of correlated orbitals (see Sec. IVI[) . KKR 
total energies are very sensitive to the angular momen- 
tum expansion used for calculation. To get accurate re- 
sults we performed LDA numerical tests up to /max = 6. 
We found that in the case of Ni and Mn to obtain con- 
verged results we need to use at least angular momentum 
expansion up to /max = 3. This cut-off was used for the 
more computationally demanding LDA-I-DMFT calcula- 
tions. 

The local problem was studied for different values of 
U in the range between 2 and 3 eV, considered accept- 
able from the results of constrained LDA calculations^ii^ 
and previous LDA-I-DMFT simulations. The tempera- 
ture was set as T = 400 K and 2048 Matsubara frequen- 
cies were used. As for the DFT part, convergence in 
the LDA-I-DMFT total energy was considered acceptable 
when the changes for subsequent iterations were smaller 
than 0.1 meV. 

At the top of Figure [3l we can see the total energy 
curves as functions of the lattice constant for the FP- 
LMTO implementation. The curves have been shifted 
with respect to their minima, so it is easier to com- 
pare them. As observed in previous calculations^, in 
DFT-LDA the equilibrium value of the lattice constant 
is slightly (3%) underestimated with respect to the exper- 
imental one. Looking at the curves for the LDA-I-DMFT 
simulations, we immediately notice that the results are 
strongly dependent on the value of the Hubbard U. 
Furthermore the best result seems to be obtained for 
U — 2 eV, i.e. for a value smaller than the widely ac- 
cepted U — 3 eV. On the other hand the curve for 
U = 3 eV seems to comprehend too strong correlation 
effects. The explanation of these results is in the pertur- 
bative nature of the SPTF solver, which tends to over- 
estimate correlation effects in fee Ni. This was noticed 
since the first implementation^, when comparison be- 
tween LDA-FDMFT resuhs with the SPTF solver and 
numerically exact Quantum Monte-Carlo solver showed 
the best agreement for U = 2 eV. Furthermore in the 
already mentioned calculation of the orbital polarization 
of Ni, it is shown that SPTF with U = 3 eV gives too 
strong correction of the orbital moment^. 

On the other hand we could be tempted to think that 
this behavior is increased by to the lack of the full self- 
consistency in the LDA-I-DMFT cycle. This doubt is 
removed by looking at the results for KKR, reported 
at the bottom of Figure [3l In fact we can barely no- 
tice any difference with respect to the energy curves of 
the FP-LMTO. It is important to empahsize how similar 
the presented results are, since the arbitrariness of the 
LDA-I-U Hamiltonian ([T]), due to the arbitrary choice of 
the correlated orbitals, is often considered as a limit of 
the orbital-dependent methods. 

Table I, where the equilibrium atomic volume Vq and 
the bulk modulus B are given, allows a more quantitative 



10 



TABLE I: Computed values of the equilibrium atomic volume Vb and the bulk modulus B for the the standard LDA-DFT 
method and for the LDA+DMFT scheme. Different strength of the local Coulomb repulsion U have been studied, at T = 400K. 
The values taken from Ref. (4^ are obtained by means of an ASA-LMTO code. 





LDA 


U = 2.0 eV 


U = 2.3 eV 


U = 3.0 eV 


GGA 


EXP 




FP-LMTO 


67.88 


76.20 


79.19 


89.48 






Vb(a.u.3) 


KKR 
Ref. 48 


66.86 
67.71 


76.28 


79.02 


85.53 


76.54 


73.52 




FP-LMTO 


260 


163 


142 


84 






B(CPa) 


KKR 
Ref. 48 


280 
270 


171 


150 


132 


186 


186 



comparison of the two implementations and with previ- 
ous DFT-LDA studies of fee Ni''^. These values of Vq and 
B have been computed with polynomial fitting of the en- 
ergy versus atomic volume curve around the minimum. 
In addition also fitting through Birch-Murnaghan equa- 
tion of state^i^ was done, leading to almost identical 
results and confirming the stability of our data. 

As for the total energy curves, the best results are ob- 
tained for U ~ 2 eV, and we can see that the inclusion 
of local correlation effects into the LDA results corrects 
both the equilibrium atomic volume and the bulk modu- 
lus in the right way. While this fact has enough interest 
by its own, we should notice that to have more precise 
results on the quantitative point of view, a more strict 
relation between solver, correlated orbitals and values of 
U is needed. Naturally it would be interesting to repeat 
those calculations with the numerically exact quantum 
Monte-Carlo solver to check if better agreement with the 
experiment can be obtained. Another interesting prop- 
erty can be deducted from the Table I: while the equi- 
librium atomic volumes are independent on the full self- 
consistency, the bulk modulus looks to be more strongly 
influenced. As expected this discrepancy is proportional 
to the strength of U. The simulation for the strongest 
value tried, i.e. U = 3 eV, shows the tendency of the FP- 
LMTO to underestimate the value of the bulk modulus 
of fee Ni. 



VIII. 7-MN 

Mn is deflnitely one of the most interesting and com- 
plex materials among pure transition metals. According 
to Hund's rule, free atom possesses a large magnetic mo- 
ment of 5/iB, and the stabilization of such large magnetic 
moments, e.g. in Heusler alloys, would represent a great 
technological advance, suitable for many applications. 

Experimentally Mn exists in four different phases. The 
low-temperature low-pressure phase is the a-phase^. It 
has a complex cubic structure with 58 atoms per unit 
cell and non-coUinear antiferromagnetic order. The lo- 
cal moment depends strongly on the atomic site, vary- 



ing between 3/is and 0, and disappears above the Neel 
temperature Tn = 95 K. At T = 1073 K there is a 
transition to the /3-phase^, a cubic structure with 20 
atoms per unit cell and small magnetic moment. Be- 
tween T = 1368 K and T = 1406 K a high-temperature 
7-phase with fee structure appears. Interestingly this 
phase can be stabilized until room temperature through 
the addition of a small amount of impurities^ or as layer- 
by-layer deposition on Cu3Au(100)^i^. Below the Neel 
temperature, about T/v = 540 K the 7-phase acquires an 
anti-ferromagnetic ground-state, which is accompanied 
by tetragonal distortion into the fct structure^i^. From 
T = 1406 K up to the melting temperature Tm = 1517 K 
there is a (5-phase, whose structure is bee and order is 
antiferromagnetic. Finally high-pressure studies have re- 
vealed a transition to an hep e-phase^ at 165 GPa. 

Such a rich phase diagram corresponds to an equiva- 
lently rich history of theoretical studies (for an extended 
and detailed review we redirect the reader to the Ref. 
[26I ). Obviously these studies have been mainly focused 
on the two "simplest" phases, 7 and 6, while the increase 
of the computational power achieved in the last ten years 
made the first ab-initio calculations of a and /? phases 
appear^i^^. 

Our main interest concerns the ground-state proper- 
ties of 7-Mn and the role of correlation effects. The de- 
scription of the electronic properties given by density- 
functional theory is undoubtedly wrong for non spin- 
polarized LDA, but it becomes more reasonable if spin- 
polarization is introduced^iiS. As for Fe, however, LSDA 
does not predict the correct crystal structure, but the 
ground-state of Mn results to be hcp^"*. Furthermore 
these strong magneto- volume effects are reflected into an 
anomalously low value of the bulk modulus'^'*. This can 
be considered as a flrst hint to strong correlation effects. 
Like for the other transition metals, the agreement of the 
calculated ground-state properties with the experimental 
data improves drastically if spin-polarized GGA is used 
as exchange-correlation potentia l^^'^^ , but the discrep- 
ancies are still the strongest of the 3d series. Further- 
more, as already pointed out by Zein^i, the anomalous 
properties of Mn do not seem to depend so strongly on 
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the magnetic phase. In fact extrapolation of experimen- 
tal data for Mn-Cu alloys to zero content of Cu shows^ 
equilibrium atomic volume and bulk modulus compara- 
ble (in a range of 10%) to pure 7-Mn, while doping by 
Cu suppresses antiferromagnetism in 7-Mn. The situa- 
tion becomes still worse if spectral properties are consid- 
ered. The only LDA-I-DMFT study available on 7-Mn 
has showni^ that inclusion of local Coulomb interactions 
is necessary for a proper description of the excitations. 
Following this work, 7-Mn seems to behave more as a 
strongly correlated metal at the metallic side of Mott 
metal-insulator transition, than as a moderately corre- 
lated metal with some deficiencies in the spectrum, as 
Ni: Hubbard bands are formed for high energies and a 
quasiparticle resonance appears around the Fermi level. 
To clarify the role of correlations and the connection be- 
tween correlations and magnetism in 7-Mn we have car- 
ried out systematic LDA-I-DMFT simulations. We have 
adopted a simple fee crystal structure in a layered an- 
tiferromagnetic phase AFMl, since previous simulations 
showed clearly this to be the equilibrium structure^s^i^. 
As already deduced in the early eighties^, the frustra- 
tion of the AFMl fee structure should imply a slight (6%) 
distortion of the lattice, but this effect has not been con- 
sidered here, since its role is not so important in compar- 
ison to local Coulomb interactions. The relation between 
correlation effects, frustration and lattice distortion will 
be the subject of future investigations. The lattice con- 
stants have been ranged from a — 6.0 a.u. and up to 
a = 7.5 a.u.. All the other computational details have 
been set as the ones used for Ni. 

The choice of the Hubbard U for Mn is not trivial at 
all, since this element was not studied as much as Ni. 
In the previous LDA-f DMFT study^- it was varied be- 
tween 3 eV and 5 eV through semi-empirical considera- 
tions. However, recent progress has been made on the 
implementation of procedures to determine the param- 
eters describing the local Coulomb interactions ah ini- 
tio. New results for the 3d transition metals have been 
obtained using the "canonical" constrained local den- 
sity approximations, and the "new" constrained random- 
phase approximation^ii^ and they locate U in the range 
2 — 4 eV for the whole series, reaching maximum val- 
ues for the half-filled systems. Given that one of these 
simulation used a basis set very similar to ours (head 
of the LMTO)^, for 7-Mn we adopted U = 2.6 eV and 
U = 3.0 eV. The corresponding Stoner parameter was 
chosen as, respectively, / — 0.8 eV and / = 0.9 cV. 

In Figure [5] the total energy curves as functions of the 
lattice constant for the FP-LMTO implementation are 
given. As for Ni, the curves have been shifted with re- 
spect to their minima to obtain a better visualization. 
From Figure 31 we immediately notice two interesting 
features in the LDA-I-DMFT total energy curves. First 
of all we can notice that, by increasing the value of U 
from zero to the accepted effective value, the minima of 
the total energy curves of the LDA-f DMFT simulations 
gradually tend to the experimental lattice constant. Fur- 
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FIG. 4: (Color on-line) Energy versus lattice constant curves 
for 7-Mn in the DFT-LDA scheme and in the LDA-I-DMFT 
scheme based on the FP-LMTO method. The zero of the 
energy of each curve is set to its own minimum value Eo 
and two chosen values of U are presented (T = 400 K). The 
lattice constant that corresponds to the experimental atomic 
volume is indicated by the arrow. In the inset we can observe 
the total energy for LDA-fDMFT simulation a.t U = 2.6 eV 
(big points) as function of the atomic volume compared to 
the standard Birch-Murnaghan equation of state (solid line). 
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FIG. 5: (Color on-line) Local magnetic moment /i and 
Galitskii-Migdal contribution to the total energy < Hu > 
as function of the lattice constant for 7-Mn. While it is 
not observable from the picture the magnetic moment of the 
LDA-I-DMFT simulation is increased with respect to its bare 
LDA value. For U = 2.6 eV the increase in the magnetic 
moment is about 0.02 /is, while for U = 3.0 eV it is about 
0.03 /is. Interestingly no magnetic moment is created if the 
starting Kohn-Sham densities is non-magnetic. 



thermore the dependency of the results from the strength 
of [/, which have been already observed for Ni, looks still 
bigger and we consider it as good indication for strong 
correlations. This impression is emphasized by another 
interesting feature noticeable from Figure ID the total 
energy curves do not appear to have a perfect parabolic 
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shape as for usual LDA or GGA simulations, or also for 
the LDA+DMFT simulations of Ni depicted in Figure 
[3l Instead they show a small kink for lattice constants 
around 6.6 a.u.. To make it more visible, in the inset of 
Figure |4] the calculated data for U = 2.6 eV are com- 
pared with a standard fitting through Birch-Murnaghan 
equation of state. This kink is a clear sign of the strongly 
correlated character of 7 Mn and reminds the one found 
in LDA+DMFT total energy curves of (5-plutonium^. In 
the latter case, there is more than just kink, there is a 
second minimum of the total energy which was associ- 
ated with the volume of monoclinic a phase. For Mn, 
there is no phase transitions with large volume jumps, 
like for Pu, but, instead, anomalies of the bulk modu- 
lus in Mn-based alloys are observedSfi.. It is important 
therefore to analyze the origin of this kink. In Figure 
[5] magnetic moments and Galitskii-Migdal contributions 
to the total energy functional are shown. We can see 
that the value of the lattice constant corresponding to 
our kink is a bit higher than the critical value for which 
the non-zero magnetic moment appears. At the onset of 
the magnetism, the competition with the local Coulomb 
interactions brings a saturation of the Galitskii-Migdal 
energy, which otherwise would be expected to decrease 
with the atomic volume (as for example we observe for 
Ni). Instead of decreasing the correlation energy, the sys- 
tem responds with an increase of the magnetic moment 
with respect to the bare LDA value. This change is so 
small that it can be barely noticed in the upper plot of 
Figure For U = 2.6 eV the increase of the magnetic 
moment is about 0.02/iB, while for U — 3.0 eV it is about 
0.03 HB- 

Given that the FP-LMTO implementation is numeri- 
cally less expensive than FP-KKR, we have made exten- 
sive calculations for 7-Mn only using the former method. 
A few simulations have been made also with the FP- 
KKR method and the same qualitative features reported 
in Figures m and O have been observed, stating again that 
for the description of the ground state properties of 3d 
transition metals the inclusion of the local correlation 
effects on the electron density is not strictly necessary. 

A more clear picture of the physical properties of 7- 
Mn can be obtained from the Table II, where equilib- 
rium atomic volume Vq, bulk modulus B and magnetic 
moment /i for our simulations have been compared to the 
exp erimental values and to the results reported in Ref. 

m 

Consistently with previous calculations, the LDA fails 
for 7-Mn and the differences with the experimental data 
are much stronger than for the other transition metals, 
e.g. Ni presented above. The atomic volume is un- 
derestimated and the bulk modulus is heavily overesti- 
mated. Moreover for 7-Mn the change of the exchange- 
correlation potential from LDA to GGA does not solve 
all the problems, and still there is an important difference 
between theory and experiments. Does the LDA+DMFT 
scheme give a better description? The simulation for the 
weakest U seems to underestimate the local Coulomb in- 



teraction. The corrections of equilibrium atomic volume, 
bulk modulus and magnetic moment are good, but they 
are too small to reproduce the experimental data. On 
the other hand the simulation for the strongest U is in 
perfect agreement with the reported values. Neverthe- 
less we must notice that the quantitative difference of the 
bulk modulus between the two LDA+DMFT simulations 
is surprisingly big. From the comparison with FP-KKR 
data, and also looking to the results for Ni, we see that 
our value is slightly underestimated because of the use 
of the basic DMFT cycle, but we can exclude that this 
effect comprehend the whole variation of B. We iden- 
tify this sensitivity of B to U as another sign of strong 
correlations. 

The reliability of the solver used in the presented cal- 
culations has been checked carefully. In fact the SPTF 
solver is a perturbative approach to the Anderson impu- 
rity model, and its application is restricted to systems 
where the Hubbard U is not bigger than the bandwidth. 
In this sense 7-Mn is a system at the border of the range 
of applicability, so that a deep investigation of the behav- 
ior of SPTF has been necessary. Given that the localiza- 
tion of the 3c? electrons depends on the atomic volumes, 
we could expect that our approximations are not valid 
for high values of the lattice constant. We surely ex- 
clude this problem since we verified that this happens 
only far away from the range of atomic volumes we were 
interested in. Another problem we could exclude was the 
fact that our approximations could simply collapse for 
all the atomic volumes driven by the strength of U. In 
fact we have studied intermediate values of U between 
U — 2.6 eV and U = 3.0 eV and all the physical prop- 
erties have shown a regular behavior, including the bulk 
modulus B. 

While we focused our analysis mainly on the anti- 
ferromagnetic phase, we tried to get more insight into the 
role of magnetism in 7-Mn through LDA+DMFT simu- 
lations of the non-magnetic phase. The results are quite 
interesting: the energy versus lattice constant curve (not 
shown here) has a regular parabolic shape with an equi- 
librium atomic volume Vq = 85.91 a.u.'^, intermediate to 
the equilibrium atomic volume of the LDA+DMFT sim- 
ulation for the antiferromagnetic phase. Obviously this is 
a consequence of the constrained zero magnetic moment, 
and no quenching of the Galitskii-Migdal energy can ap- 
pear. The increasing strength of the correlation energy 
is observable also in a huge drop of the bulk modulus 
with respect to its bare LDA value: B — 57 GPa, per- 
fectly consistent with the already mentioned experimen- 
tal data for 7-MnCu alloys^-, after extrapolation to zero 
content of Cu at room temperature. As before we have 
checked whether the SPTF solver is applicable or not 
to our system. We have found that our approximations 
lose validity for atomic volumes larger than 100 a.u.'^: 
the localization effects are heavily overestimated and the 
crystal tends to collapse into an atomistic system. Fortu- 
nately this threshold is well above the equilibrium values, 
so that we can still consider our results as reliable. 
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TABLE II: Computed values of the equilibrium atomic volume Vb, the bulk modulus B and magnetic moment /i of 7-Mn for 
the standard LDA-DFT method and for the LDA+DMFT scheme. Different strengths of the local Coulomb repulsion U have 
been studied, at T = 400K. The values taken from Ref. [2^ are obtained by means of a USPP-PAW (ultrasoft pseudopotential 
projector augmented plane-wave) code, and using the Murnaghan equation of state^i^. The experimental values for the 
atomic volume and the magnetic moment come from Refs. [56ll69l . and are obtained as extrapolation to room temperature of 
high temperature data. The values of the bulk modulus are more uncertain and come from Refs. [2^l70l . 





LDA 


U = 2.6 eV 


U = 3.0 eV 


GGA 


EXP 


Voia.u.^) 


FP-LMTO 
Ref. 25 


69.18 
68.36 


81.17 


88.61 


82.32 


87.30 ^ 87.60 


B(GPa) 


FP-LMTO 
Ref. 25 


313 
310 


213 


88 


95 


90 -^ 130 




FP-LMTO 
Ref. 25 


0.00 
0.00 


1.74 


2.30 


2.40 


2.30 



IX. CONCLUSIONS 

In this paper we have presented two different total en- 
ergy implementations for the LDA-I-DMFT method, us- 
ing the SPTF solver for the solution of the local prob- 
lem. Our codes have been tested through the study of the 
ground-state properties of fee Ni. The results have been 
very encouraging, showing good agreement with experi- 
mental data and in a particular a weak dependence on 
the implementation or on the choice of the local orbitals. 
Furthermore a tendency of the basic LDA-I-DMFT cycle 
to underestimate the bulk modulus with respect to the 
fully self-consistent cycle has been observed. 

The main scientific aim of this paper has been the 
analysis of the role of local correlations in 7-Mn. Clear 
signs of strong correlations have been found and the 
LDA-fDMFT method has been shown the ability to 
treat the non-magnetic and anti-ferromagnetic phases on 
the same footing, improving considerably the results ob- 
tained with usual one-particle approximations. 

Finally the results presented here stimulate future re- 
search. The main question concerns the origin of a kink 
in the total energy curves and the role of the tetrago- 
nal distortion of the fee lattice on the correlation effects 



of the antiferromagnetic phase of 7-Mn. This last study 
can be particularly interesting for the calculation of the 
elastic properties. In addition the influence of the choice 
of the solver on the description of 7-Mn needs more in- 
vestigation. 
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